A search for low-mass WIMPs with EDELWEISS-II heat-and-ionization detectors 
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We report on a search for low-energy (E < 20 keV) WIMP-induced nuclear recoils using data 
collected in 2009 — 2010 by EDELWEISS from four germanium detectors equipped with thermal 
sensors and an electrode design (ID) which allows to efficiently reject several sources of background. 
The data indicate no evidence for an exponential distribution of low-energy nuclear recoils that 
could be attributed to WIMP elastic scattering after an exposure of 113 kg-d. For WIMPs of mass 
10 GeV, the observation of one event in the WIMP search region results in a 90% CL limit of 
1.0 x 10~ 5 pb on the spin-independent WIMP-nucleon scattering cross-section, which constrains 
the parameter space associated with the findings reported by the CoGeNT, DAMA and CRESST 
experiments. 

PACS numbers: 95.35.+d, 14.80.Ly, 29.40.Wk, 98.80.Es 
Keywords: Dark Matter, Cryogenic Ge detectors, WIMP searches. 



I. INTRODUCTION 

Within the current cosmological concordance model, 
a large fraction of the mass content of the observ- 
able universe is made of dark matter, as demonstrated 
by a range of observations but the nature of which is 
still unknown |1]. In particular, the dynamics of our 
galaxy can be explained by the presence of a dark mat- 
ter halo with a solar neighborhood density of 0.3 ± 
0.1 GeV/cm 3 . Weakly Interacting Massive Particles 
(WIMPs) constitute a generic class of dark matter can- 
didates. The calculated thermal relic density of WIMPs 
roughly matches the measured density on cosmological 
scales. Although several constrained WIMP models point 
to WIMP masses M x ~ 100 GeV, lighter fermions with 
masses down to ~ 2 GeV are possible [113. 

The direct detection of low-mass (~10 GeV or lower) 
WIMPs is challenging because the recoil energies gener- 
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ated by the elastic scattering on nuclei of such low-mass 
particles are close to the experimental thresholds of ex- 
isting detectors. Experiments using a range of technolo- 
gies have either constrained 0, [f| or found hints [fUl] 
of such WIMPs with cross sections for spin- independent 
scattering on a nucleon ctsi of the order of magnitude of 
10~ 4 pb. There is a discussion regarding these results due 
to uncertainties on efficiencies, energy scales and residual 
backgrounds at low recoil energies (see for example Q). 

In this article we present a search for low-mass WIMPs 
using data collected by the EDELWEISS-II detectors 
in 2009 - 2010. A WIMP search optimized for WIMP 
masses above 50 GeV has already been published [10l - fl2j , 
where data from ten heat-and-ionization germanium de- 
tectors with interleaved electrodes were used. The analy- 
sis threshold was set at a recoil energy of 20 keV in order 
to ensure a maximum exposure in a recoil energy range 
where the behavior of all of the ten detectors was homo- 
geneous and well understood, both in terms of expected 
background and rejection capabilities. This strategy, 
however, is not adequate for a search dedicated to models 
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in which the WIMP mass is of the order of 10 GeV, for 
which the highest expected recoil energy is of the order of 
10 keV. In the study presented here, we use a restricted 
data set, selected on the basis of detector thresholds and 
backgrounds, for which a low-background sensitivity to 
nuclear recoils down to 5 keV could be achieved. By lim- 
iting this analysis to energies less than 20 keV, the results 
presented in this article are completely independent from 
the previously published limits [ll( . After a presentation 
of data selection, we discuss the potential residual back- 
grounds in the WIMP search region. A total exposure 
of 113 kg-d is obtained, from which we derive constraints 
on ctsi for WIMP masses in the range 7 — 30 GeV. 



II. ANALYSIS STRATEGY 

We use data collected during 14 months in 2009 — 2010 
by an array of ten ID detectors within the experimental 
setup described in detail in Ref. [Ill- In each bolometer, 
individual interactions are detected by a thermal (heat) 
sensor, and the charge created during an interaction is 
collected by a set of six electrodes - two fiducial elec- 
trodes, two veto and two guard channels. The electric 
field geometry created by the voltages applied to the 
electrodes isolates a central fiducial volume within the 
detector. Interactions within the fiducial volume pro- 
duce equal amplitude and opposite polarity signals on 
the fiducial electrodes, while surface interactions induce 
signals also on the veto and guard electrodes [l3j . 

To be sensitive to recoil energies well below 20 keV and 
to minimize backgrounds in the WIMP search region, we 
rejected, a priori, six detectors among the ten bolometers 
available. Four of these six detectors, which had non- 
functioning electrodes and/or poor resolutions for one or 
several channels, were rejected to guarantee the discrim- 
ination of the gamma-ray and surface event backgrounds 
around 10 keV. Another detector was rejected due to the 
presence of a relatively intense accidental source of 210 Pb 
in its vicinity, leading to a large number of low-energy 
and very low ionization-yield events induced by a radi- 
ation and Pb recoils. Finally, one detector was rejected 
because a low-energy gamma-ray background within its 
fiducial volume, which was 4 times larger than in other 
bolometers was observed El- 
Periods of good data quality were selected in order to 
provide a homogeneous dataset for each of the four de- 
tectors used in the analysis. This time period selection 
was based on the measured baseline resolutions of each 
channel on an hour-by-hour basis. A large part of the 
corresponding exposure loss is due to noisy periods in 
the heat channels. That noise can be associated with 
certain readout channels that exhibit stronger than ideal 



coupling to the microphonic noise caused by the cryo- 
genic equipment that was in operation during data tak- 
ing. Depending on the detector, the number of live days 
selected by these cuts ranges from 120 to 247 days. Based 
on the achieved thresholds on heat and ionization signals 
described below, it is expected that a significant fraction 
of the sensitivity to low-mass WIMPs will be driven by 
one of the detectors, called ID3, for which 197 live days 
were accumulated with average FWHM resolutions on 
the heat and fiducial ionization of 0.8 keV and 0.7 keV 
(electron recoil energy scale), respectively. The measured 
fiducial mass of these 400 g ID detectors is 160 g [To| . 
This measurement was done using single-scatter events 
from cosmogenic activation gamma-ray lines with ener- 
gies of 9.0 and 10.4 keV. The fiducial volume is primarily 
related to the electric field configuration inside the crys- 
tals, and the fiducial cuts are based on signals from veto 
electrodes, which do not depend on the fiducial energy. 
As a consequence, the measurement is applicable in the 
energy range considered for this study, and yields a total 
WIMP search exposure of 113 kg-d. 

The event reconstruction is identical to that carried 
out in [llj], as are standard event-based quality cuts. In 
particular, a cut is applied on the x 2 of the heat pulse 
reconstruction. The efficiency of this cut is independent 
of heat energy, and is 98.7%, determined from a mea- 
surement using the aforementioned internal radioactiv- 
ity lines at 10 keV. For the WIMP search we also re- 
ject events in coincidence with interactions recorded by 
any other surrounding bolometer in operation, or in the 
muon veto. The deadtime associated with this rejection 
is 0.2 %. Finally, we restrict the WIMP search to mea- 
sured recoil energies below 20 keV, making this data set 
independent from the one analyzed in |ll| . 

In order to search for nuclear recoil (NR) energy signals 
near the threshold, it is advantageous to evaluate the 
recoil energy of each event using an estimator with a very 
good energy resolution. We evaluate the recoil energy E r 
under the assumption that the event is due to a nuclear 
recoil as in [4J by inverting the following formula: 



l + \V\/3 
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—Q n (E r ) 



1 However, it has been checked that including data from this de- 
tector does not significantly affect the resulting combined WIMP 
limit. 



This describes the relationship between the heat energy 
E)x, calibrated by gamma-ray interactions, and the recoil 
energy E r , taking into account the heating of the detector 
associated with the charge drift under the polarization 
voltage V in volts (Luke-Neganov effect [H], [ll|). Here 
Qn{E r ) = 0.16 (EV/keV) ' is the nuclear recoil ioniza- 
tion quenching factor used in [llij . and validated down 
to 5 keV using neutron calibrations. This estimator pro- 
vides a better energy resolution on E r for nuclear recoils 
than the estimator used in [ll|, [TtJ ■ Note that for a fidu- 
cial electron recoil, E r overestimates the real energy by 
a factor ~ 2 for the typically applied voltages. 

At low energy an efficiency loss appears due to the on- 
line trigger. Events are recorded based on upward flue- 
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tuations of the bolometer temperature measurement (fil- 
tered to remove noise). The trigger threshold was varied 
during data acquisition, depending on the noise condi- 
tions of the channel, and was recorded to disk. At any 
given time, the dependence of the trigger efficiency on 
E r can be described as a function of recoil energy as 
^online {E r ) = 0.5[f + Erf((£ r - E thlcsh )/aV2)}, where 
-Ethresh and a are the recorded threshold and the mea- 
sured resolution, in NR energy scale. This parametriza- 
tion was validated using the flat, low-energy Compton 
plateau in gamma-ray calibration data. For the selected 
data, the average trigger efficiency is 78% at 5.0 keV (NR 
scale) and 90% at 6.3 keV. 

We now consider signals from the ionization channels, 
which are exclusively used to reject the main backgrounds 
of nonfiducial interactions such as surface beta radioac- 
tivity, gamma-ray-induced interactions in the fiducial 
volume, and ionizationless events. 

• Fiducial events are selected by requiring the ab- 
sence of any signal on the veto and guard elec- 
trodes. Since a significant fraction of the ioniza- 
tion signal is always present on these channels for 
surface events, this cut efficiently rejects nonfidu- 
cial interactions down to low energy. Although the 
rejection performance depends on the amount of 
energy detected in the nonfiducial electrodes, the 
efficiency loss induced by the cut does not depend 
on energy. The effect of this loss is included in the 
measured fiducial mass [Toj . 

• The reconstructed nuclear recoil energy E r , and 
the weighted fiducial ionization energy Ei, in elec- 
tron recoil scale, are calculated for events that 
pass the selection described above. Ei is the com- 
bination of fiducial electrode signals which results 
in the best ionization resolution. In the (E r ,Ei) 
plane, the fiducial gamma-ray interactions are lo- 
cated in a band centered along the line Ei = 
E r (l + Q n |V|/3)/(l + |V|/3), the width of which is 
determined by the heat and ionization resolutions. 
This background is rejected by selecting events be- 
low the band in which 95% of the gamma-rays are 
located. 

• In WIMP search conditions, a large number of 
events were recorded with signals only on the heat 
channel. These ionizationless events have several 
sources: Pb recoils from surface alpha radioactiv- 
ity, internal radioactivity of the NTD sensors, and 
non-Gaussian noise of electromagnetic and thermal 
origins. In order to reject efficiently these events, 
we take advantage of the redundancy of both fidu- 
cial electrode measurements and require the pulse 
fit algorithm to reconstruct signals from these two 
electrodes within 30 ^s from each other in the 
recorded traces. In addition, only events with a 
weighted fiducial ionization energy Ei larger than 
twice the corresponding FWHM are selected, the 
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FIG. 1. Top panel: Distribution of the fiducial ionization en- 
ergy, Ei, in neutron calibration for the ID3 detector, before 
(black) and after (red) the ionization cut. The spectrum con- 
tains mostly nuclear recoils and a tail of ionizationless events. 
Bottom: cut efficiency determination. The fit function (con- 
tinuous red line) is e ion = 0.95 [1 - exp(-1.87(£ ion - 1-25))] 
where Ei on is in keVee (electron recoil scale). The dashed 
lines represent the efficiency functions estimated in the same 
way for the three other detectors used in this search. 



average value of which ranged from 1.4 to 1.9 keV, 
depending on the detector. 

The dependence of this last cut on the reconstructed 
fiducial ionization is measured from neutron calibrations, 
which provide a large sample of low-energy, WIMP-likc 
pulses with ionization noise conditions identical to WIMP 
search conditions. Figure Q] presents the distribution of 
Ei in neutron calibration for the ID3 detector, before and 
after this ionization cut. The measured efficiency func- 
tion £ion(Ei) is obtained through an analytic fit to the 
ratio of both histograms. Because of the presence of a 
residual tail of ionizationless events even in the short run 
time with a neutron source, the measured efficiency is 
underestimated at low energy, and therefore provides a 
conservative estimate of the WIMP sensitivity. 

The WIMP signal density is calculated based on the 
WIMP-induced nuclear recoil spectrum as a function of 
true recoil energy po(E r o) as parametrized in [18| . For 
a given detector and a given WIMP mass, taking into 
account the measured resolutions of fiducial ionization 
energy <jj (electron recoil scale) and heat energy (NR 
scale) oy, the WIMP signal density in the (E r ,Ei) plane 
is then: 

p(E ^ Ei) = 6onlinc(gr) £i on(gi) f rf^^^) 

lna r CFi J 
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(E r — E r0 ) 2 (Ei — Q n E r0 ) 2 



This function is normalized for a cross section ergi = 
10~ 6 pb and an exposure of 1 kg-d, as it scales triv- 
ially with these parameters. For each detector and each 
WIMP mass, we then define a WIMP search region in 
the (E r , Ei) plane as the region containing 90% of all the 
calculated WIMP signal density below the gamma rejec- 
tion cut, mentioned above. For example, Fig. |3] shows 
the function p and the WIMP search region for ID3 and 
M x = 10 GeV. 

The parameters Q n (E r ), a r and <7j, entering in the 
definition of this WIMP search region, were cross-checked 
using calibrations down to the lowest relevant energies: 
the width of the gamma-ray band as well as the position 
and width of the neutron band are compatible with the 
measured baselines and with the parametrization Q n — 
0.16 (£ r /keV)°- 18 from 0. 



III. BACKGROUNDS, RESULTS AND 
DISCUSSION 

Several backgrounds were anticipated : 
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FIG. 2. Scatter plot of heat vs ionization for all selected events 
in the 113 kg-d WIMP search. Neutron calibration data are 
also displayed as grey dots. The continuous (dashed) blue line 
represents the 95% gamma-ray rejection cut for detector ID3 
(ID401). The blue, circled points are events contained in the 
WIMP search region for at least some WIMP masses. 



• Ionizationless pulses were the most prominent 
background at very low energy, since we observe 
between 1000 and 5000 such events above 5 keV 
per detector. The rejection factor of the ionization 
cut is of the order of ~ 10 -6 , resulting in a negligi- 
ble contribution to the remaining background. 

• The residual fiducial gamma-ray background is es- 
timated by extrapolation of the observed rate of 
gamma-ray events in the energy range for which the 
WIMP search region is limited by the 95% gamma 
rejection cut. This range depends on the detector 
and WIMP mass under consideration. For exam- 
ple, for M x = 10 GcV its average value is 5.9 — 
8.8 keV. The estimated background from Gaussian 
leakage is then 2.5% times the measured number 
of gamma events in that energy range, which leads 
to an overall estimate of (1.2 ± 0.2) events. For 
M x = 30 GeV this figure becomes (1.1±0.2) events. 

• An upper limit on neutron backgrounds from dif- 
ferent sources, as determined by the same combi- 
nation of simulations and measurements as in [IT| , 
was found to be 1.7 events in the energy range 
5 — 20 keV for the exposure of 113 kg-d. A most 
probable contribution of 1.0 events is expected from 
radioactivity in the warm electronics, cables and 
connectors. 

• Surface interactions due to beta radioactivity are 
rejected but we have not measured rejection factors 
in the relevant 5—15 keV energy range. 
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FIG. 3. Color map of the WIMP signal density p(E r , Ei) to- 
gether with the 90% WIMP search box (red line) for a 10 GeV 
WIMP in detector ID3. Overlaid are background (black dots) 
and neutron calibration data (grey dots) for the same detec- 
tor. The average position of the fiducial gamma-ray band and 
the corresponding 95% rejection cut are shown as continuous 
and dashed blue lines, respectively. 
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FIG. 4. 90% CL Poisson limit on <t S i as a function of WIMP 
mass derived from the analysis of the four bolometers (bold 
red line). We also show the location of potential WIMP sig- 
nals from the CoGeNT @], CRESST @] and DAM A [H[ ex- 
periments, as well as constraints from EDEIWEISS-II [llj, 
CDMS-II [13] and XENON100 [U, and the dedicated low- 
mass searches by CDMS-II 0] and XENON10 0. 



Figure [2] shows the (E r ,Ei) scatterplot of the events 
recorded in the 113 kg-d exposure selected for the WIMP 
search. Most events are compatible with gamma-rays in- 
teracting in the fiducial volume. The distribution of Ei 
for these events shows the presence of the cosmogenic ac- 
tivation lines at 9.0 and 10.4 keV, as well as less intense 
lines in the 5 — 7 keV energy range. While a few events 
are present in the nuclear recoil region, there is no indi- 
cation for an exponential distribution of nuclear recoils, 
in particular at energies below 10 keV, where most of the 
WIMP signal for M x ~ 10 GeV is expected. Figure 
shows the signal density p(E r ,Ei) and the WIMP search 
region corresponding to M x = 10 GeV for the detector 
ID3, showing the absence of a nuclear recoil signal in 
spite of the sensitivity to such a signal demonstrated by 
neutron calibration data. 

The total residual background in WIMP search regions 
ranges from one event for M x = 10 GeV to three events 
at 30 GeV. This is compatible with the expected back- 
grounds described above. For relatively large masses, the 
sensitivity of this study is therefore background limited, 
while for low WIMP mass we are primarily limited by 



the achievable threshold. 

To derive the corresponding limit on <tsi as a function 
of WIMP mass, we count the events located within the 
WIMP search region defined previously, simply adding 
the number of events over the four selected detectors. 
The expected number of events for a given value of 
<7si is obtained by integrating the WIMP signal density 
p(E r ,Ei) over the WIMP search region, and summing 
over the detectors. We then compute a 90% CL limit 
on the total number of events using Poisson statistics. 
The obtained limit is shown as a function of M x for 
7 <M X < 30 GeV in Fig. g] 

This result extends the sensitivity of the previous 
EDELWEISS-II analysis down to WIMP masses below 
20 GeV. We exclude the entire zones corresponding to 
the interpretation of DAMA @ and CRESST Q results 
in terms of elastic, spin- independent WIMP scattering. 
These zones extend down to a WIMP mass of 9 GeV, 
which is the "WIMP-safe" mass [19[ where our experi- 
mental sensitivity is around 1 % of the total WIMP sig- 
nal recoil spectrum. While we significantly constrain part 
of the CoGeNT parameter space Q, we cannot exclude 
the region corresponding to M x < 8 GeV due to lack of 
sensitivity to nuclear recoil energies below 5 keV. 

Our results provide a cross-check of those obtained 
with other technologies, using duahphased xenon TPCs 
(@)H3j) or germanium detectors ([j, In contrast to 
the latter, the WIMP search with EDELWEISS ID detec- 
tors presented in here is not limited by the presence of a 
large background as seen in CDMS and CoGeNT, which 
may be due to surface interactions. Although we do not 
have a precise measurement of the rejection factor for this 
background below 15 keV, the WIMP search data pre- 
sented here demonstrates that it remains excellent down 
to 5 keV. This provides strong motivation for measuring 
this rejection factor with new-generation detectors. We 
aim to make further improvements in lowering energy 
thresholds and increasing the mass while maintaining a 
background-free region of interest for low-mass WIMP 
searches, including annual modulation measurements. 
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